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We experimentally confirm the functionality of a coupling element for flux-based superconducting 
qubits, with a coupling strength J whose sign and magnitude can be tuned in situ. To measure 
the effective J , the groundstate of a coupled two-qubit system has been mapped as a function of 
the local magnetic fields applied to each qubit. The state of the system is determined by directly 
reading out the individual qubits while tunneling is suppressed. These measurements demonstrate 
that J can be tuned from antiferromagnetic through zero to ferromagnetic. 
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One approach to building useful quantum processors 
is to abandon the gate model in favor of physics-inspired 
models, such as adiabatic quantum computation (AQC) 
and its variants [1]. In AQC one adiabatically evolves 
the processor from the groundstate of an initial Hamilto- 
nian Tii (chosen such that this groundstate can be read- 
ily prepared) to the groundstate of a final Tif, encoding 
the solution to the problem of interest. For example, an 
Ising magnet is a physical system that could be harnessed 
for this purpose [2]. AQC is known to be computation- 
ally equivalent to the circuit model of quantum computa- 
tion [3]. Realistic AQC architectures exist [4] that could 
be used to produce accurate approximate solutions to 
NP-complete problems [5]. One such architecture, based 
on superconducting electronics [6], requires devices that 
couple qubit pairs with in situ programmable coupling 
magnitude and sign [7]. In this letter we demonstrate 
such a coupler between two superconducting flux qubits. 

The qubit design used is a bistable rf-SQUID, magnet- 
ically biased near its degeneracy point by an external flux 
~ ('i'o = h/2e is the flux quantum) [8, 9]. The 
qubits are denoted by a, c in Fig. 1, with bias controls 
/°, /° (/* = $!p/$o)- Compound Josephson junction 
(CJJ) loops d, e with biases are employed to tune 

the critical currents [10]. For readout, the qubits are 
inductively coupled to their own dedicated dc-SQUIDs 
/, g [11], biased to points of high flux sensitivity via a 
shared control //'^. The switching currents were mea- 
sured by ramping the bias currents if, ig and monitor- 
ing the voltages at points Vf, Vg. Tunable intcrqubit 
coupling is mediated via a monostablc rf-SQUID h with 
separate flux bias [12]. This approach and that of 
Refs. [13] and [14] are recent demonstrations of in situ 
sign- and magnitude-tunable coupling elements between 
flux qubits. 

We couple the qubits to the coupler loop via trans- 
formers leading to mutual inductances M°*, M'"^ so that 
the qubits are influenced the coupler's persistent current 
II = 4^sin(27r/*'), not the control flux Here /^ the 
total flux in loop h, is a function of [15]. The persis- 
tent current 7° in qubit a alters the flux applied to h by 



jaj^ab ^ rpj^jg additional flux alters Ip by an amount 
SI^ w /^M°*'x^ where the susceptibility x'' = dl^/d^l. 
Sip in turn alters the flux applied to qubit c by SI^M^'^, 
causing an interaction of energy 



J = M'"'M'"'x Iplj 



(1) 



where 1^ is the persistent current in qubit c. Thus 
the coupler mediates an effective qubit-qubit interac- 
tion. Since 1^ and therefore x'' S'^e periodic in with 
alternating sign, this interaction can be ferromagentic 
(FM,x'' < 0), antiferromagnetic (AFM,x'' > 0), or zero. 

The energy spacing between the ground and first ex- 
cited states of the coupler was designed to be much 
higher (^ 10 GHz) than typical qubit splittings. As 
such the coupler is expected to remain in its groundstate 
and the effective low-energy qubit-coupler-qubit Hamil- 
tonian [12] can be written as follows: 



J2 (6V('''+A%(«))+J(/, 



(2) 



q—a,c 



where and ci'^ are Pauli matrices for qubit q. Here 
e"^ and A'' represent the energy bias and tunnel splitting 
for the individual qubits, as indicated in Fig. Ic. Note 
that cx /|/« and J [Eq. (1)] scale with I^ and In 
turn, I^ and A"? are functions of and /| for q = a 
and c, respectively. Tuning these latter two flux biases 
provides a means of annealing the system from a quan- 
tum regime, where there is appreciable tunneling, to the 
classical regime where A'' 0. 

The circuit was fabricated on an oxidized Si wafer us- 
ing a Nb trilaycr process with wiring layers isolated by a 
sputtered Si02 dielectric [16]. The qubit inductance was 
designed to be 500 pH. The total parallel capacitance of 
each CJJ (single junction area A = 0.7 x 0.7 fim^) was 
designed to be C = 20 fF. Fitting of macroscopic res- 
onant tunneling peaks [17] for a similar device at large 
/9 yielded C« = (33 ± 3) fF. The discrepancy may be 
due to capacitive loading of the junctions by nearby 
wiring. The unsuppressed /| of each CJJ was designed 
to be 2.5 /iA. The single-junction {A = 0.6 x 0.6 ^m^) 
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FIG. 1: a) Device schematic: SQUID loops i, mutual induc- 
tances M'-', flux bias lines current bias lines ii and volt- 
age taps Vi as indicated (see text), b) Photograph of device, 
c) Energetics of the coupled qubit system. Qubits are repre- 
sented as bistable potentials with pseudospins as indicated. 
Local biases e' are applied to each qubit. The lowest lying 
qubit states are separated by SE^ = ^(ei)^ + (A»)2. The 
coupler is represented as a monostable potential that medi- 
ates the interaction J [see Eq. (2)]. 



coupler has designed inductance L' 

Tb — no,, A rr^rinrr «b = O-T^rfcrfc 



: 240 pH and 

= 0.9 /iA, giving 13" = 27ri''4V$o = 0.65. Fits to the 
measured M'^'^Ip and M'"^Ip (explained below) yielded 
= 0.63± 0.02. The mutual inductances were designed 
to be M"'' = M''^ = M"?^ = 25 pH. Measurements on 
similar breakout structures yielded M'^^ = (28 ± 2) pH. 

Qubit states were inferred from the bias currents if, 
ig required to switch the dc-SQUIDs into the voltage 
state. These elements had designed maximum critical 
currents of 2.5 /xA for //'^ = (A = 0.7 x 0.7 /xm^). The 
bias currents were ramped linearly from zero to 2.75 /iA 
in 80 fjs (see Fig. 2) and voltages monitored at points 
Vf and Vg. Voltage trigger thresholds were set at ~ 1.5 mV 
and monitored with 10 ns timing resolution. The two per- 
sistent current states of each qubit generated a flux dif- 
ference of (16 ± 1) m$o in their respective detector with 
the width of the switching distributions limited by the 
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FIG. 2: Bias sequence. CJJ biases /; 
held for 60 fis, and then returned to zero over 20 /iS. Readout 
dc-SQUID flux bias is held at zero during qubit evolution 
and raised to its operating point during the current bias ramps 
on if, ig. Control lines and are held constant 

throughout this procedure. 



detector sensitivity of 1.2 m$o- Using the design value 
M^f = IVF-a = 11 pH, we estimate 2/9 = (3.0 ± 0.2) ^i A 
when = = 0. 

The device was mounted in an Al box in a dilution 
refrigerator with base temperature below 10 mK. Two 
coaxial cryoperm shields surrounded the sample area. 
Battery powered current control electronics were located 
in an rf shielded room together with the fridge, and inter- 
faced though fiber optics with a PC outside. All control 
lines had discrete element filters at the 1 K stage and mix- 
ing chamber as well as copper powder filters attached to 
the latter. All flux bias couplings were weak, minimiz- 
ing the dissipation introduced into the qubits from the 
environment via the wiring in the fridge. 

In order to account for cross coupling between the flux 
bias controls we measured the mutual inductance matrix 
between lines in situ by tracking features with known 
period $o in each of the qubits and the coupler. Better 
design can avoid linear crosstalk but not nonlinear qubit 
biasing due to 1^ [12]. To calibrate the latter we sup- 
pressed of one of the qubits q and measured the state 
of qubit q' ^ q versus and /J . The degenerate bias 



$9 now follows from M't'^I^ + ^ i$o- Fitting these 
data determined /3'' as noted above. In addition such 
measurements were used to generate Af ^ ''x^ which can 
then be used to predict the functional form of J (Eq. 1) 
up to a prefactor Mi^I^I^. 

We measured J by mapping the interacting ground- 
state versus e°, e'^ in the limit A' | J| where Eq. (2) 
is nearly a classical Ising Hamiltonian with eigenstates 
|ac) = ITT), iTi), UT), and |ii). The groundstate is 
found via an annealing procedure: are first increased 
by biasing the CJJs at f^''^ = ^ to render the qubits 
monostable, thus initializing them in a known ground- 
state. Thereafter A'^ are slowly lowered back to their 
minima at /^'"^ = (see Fig. 2). The system will re- 
main in the groundstate during the lowering of A'^ un- 
til the rate of evolution of H. exceeds a limit set by the 
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FIG. 3: Measured probabilities for each of the four possible 
flux states for small A,, at coupler flux bias = —0.03. 



physics of a Landau-Zcncr (LZ) transition [18]. The mea- 
sured fluxes represent the groundstate of Eq. (2) at the 
last instance before this transition. Thus according to 
Eq. (1), the observed J will depend on particular values 
of Ip{f^) and Ipif!^:), to be determined experimentally, 
which can be less than the maximum values observed at 
- /I = 0. 

Coupler performance was demonstrated at three values 
of corresponding to maximum AFM (/^ = 0), zero 
[/' = /j=o G (0, i)] and maximum FM (/,^ = i) cou- 
pling. At each bias point, the above annealing sequence 
was repeated 1024 times and the probabilities J'(ltT)) 
etc. of finding the four classical flux states were deter- 
mined. For this, direct readout of individual qubits in 
our design is a key advantage. As an example, Fig. 3 dis- 
plays P{\ac)) for = —0.03. The flux states observed 
near the corners (/^ — \,fx~\) = (±4, ±4) x 10~^ agree 
with those expected from Eq. (2) for e°, e'^ » J. While 
P(liT)) (Fig. 3a) and P{\ii)) (Fig. 3c) exceed 95% at the 
top and bottom left corners respectively, it appears that 
P(ITT)) (Fig. 3b) and P(|Ti)) (Fig. 3d) only reach ~85%. 
Further investigation revealed that this was due to mea- 
surement crosstalk from the readout of qubit c. Note 
that qubit c is read first in the bias sequence depicted in 
Fig. 2. It was verified that qubit a could be read without 
interference if the ig control pulse were absent. However, 
switching of the readout SQUID g into the voltage state 



was found to disturb the state of qubit c, which in turn 
influenced the state of qubit a prior to its measurement. 
Observations indicate that the |ti) and |tt) states were 
most susceptible to corruption on this particular device. 

Groundstate stability diagrams were generated by de- 
termining which flux state occurs with the highest prob- 
ability at each point (/°,/^)- The experimentally de- 
termined boundaries are shown in black in Fig. 4. In 
Fig. 4a, fl = -0.03, and the AFM states |Ti) and |it) 
dominate the groundstate map. One clearly observes a 
boundary between these two states, which only occurs if 
,/ > 0. This demonstrates that our device can provide 
an AFM qubit-qubit interaction. Next, for = 0.30, 
Fig. 4b shows the critical case of zero coupling (J = 0): 
the state of each qubit is independent of the flux applied 
to the other, resulting in a cross structure with the four- 
fold degeneracy point at (/°,/^) = (^,5)- Finally, for 
0.52 (Fig. 4c), the FM regions |U) and |tT) are en- 
hanced. The new boundary between them is present only 
if J < 0, indicating a FM qubit-qubit interaction. Ide- 
ally in all cases the boundaries should have intersected 
the J = degeneracy point. However, the transitions 
appear to have translated by {fx,fx) (0,— 0.1)m<I>o 
in the AFM case and by (/^, /^) ~ (-0.5, -1.5) m<I>o m 
the FM case. Further investigation revealed that this dis- 
crepancy was due to the readout crosstalk problem noted 
earlier and was aggravated by strong interqubit coupling 
in the FM case, which then altered the population statis- 
tics. 

We determine J from either the displacement be- 
tween the |iT):|TT) and |ii):|Ti) boundaries (vertical 
dashed fines in Fig. 4) or between the |tT):|Ti) and 
liT):Ui) boundaries (horizontal dashed fines). According 
to Eq. (2), in the limit A'/IJ] < 1, the former should 
yield a net flux equal to 2\J\/Ip and the latter 2| J|//p. 
The observed horizontal and vertical displacements are 
equal to within experimental error for both the FM [mean 
(0.6 ± 0.2) X 10-3] and AFM [mean (2.0 ± 0.4) x 10-^] 
cases. One concludes that the qubits and transformers 
are reasonably symmetric, i.e., J" w = 

The measured coupled flux values J(/^)/(/^<I'o) 
are shown in Fig. 5 together with the measured 
diM'i^ II) / d^l scaled by M'?''/^, where is used as a 
free parameter. The observed amount of coupled flux is 
in reasonable agreement with the predicted form (Eq. [1]) 
if the qubit currents are ~ = (0.5 ±0.1) /xA at the last 
moment before the LZ transition. Decreasing the mag- 
nitude of the ramp rate of and /| should allow the 
system to evolve more slowly as A' decreases, thus shift- 
ing the inevitable LZ transition to a larger I^. As such, 
this measure of /| at which adiabatic evolution appears 
to terminate is not due to a fundamental limitation. 

The key result of this study is a clear demonstration of 
a sign and magnitude tunable coupler between two super- 
conducting flux qubits. While originally designed for a 
specific quantum processor architecture, this device may 
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FIG. 4; Groundstate stability diagrams: a) AFM coupling, b) zero coupling and c) FM coupling. The calibrated degeneracy 
lines for the individual qubits and extreme values of coupled flux are indicated by solid and dashed lines, respectively. 
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FIG. 5: Coupling strength (in units of coupled flux) versus /j! 
as obtained from groundstate mapping and from dM''''I^/d$l 
scaled by M"*/' with = (0.5 ± 0.1) M [see Eq. (1)]. 



prove useful for various superconducting electronics ap- 
plications, both quantum and classical, requiring in situ 
tunability. Further, this work demonstrates the utility of 
a multi-qubit readout technique, directly measuring indi- 
vidual qubits, which is scalable to larger qubit numbers 
and may be of importance in future practical quantum 
processors. 
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